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Abstract 

QCD sum rules using polynomial kernels are used to evaluate the 
strong part of the proton-neutron mass difference 5M np in a model in- 
dependent fashion. The result for the mass difference turns out to depend 
sensitively on the value of the four quark condensate ((qq) 2 ) and repro- 
duces the experimental value of 8M np for ((gg) 2 ) ~ 2{qq} 2 . 
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1 Introduction 



The QCD sum rule method introduced by Shifman, Vainshtein and Zakharov |T] 
has extended the applicability of QCD far beyond simple perturbation theory. 
The method was adapted to the case of nucleons by Ioffe [2] and independently 
by Chung, Dosch, Kremer and Schall [3J. These authors showed how to approach 
one of the fundamental problems of hadronic physics, the calculation of the 
baryon masses from the Lagrangian and the vacuum condensates of QCD. 

Several authors undertook the study of the neutron-proton mass Sm np difference 
using the external field method of Ioffe and Smilga[2j .These authors use Borel 
(Laplace) transforms of correlators in order to suppress the unknown contribu- 
tions of the continua which appear in dispersion integrals. 

I present here a different approach which uses instead as kernels in the dis- 
persion integrals simple polynomials the coefficients of which are determined 
by the general shape of the physical spectrum and which, as will be shown, offer 
stronger damping of the unknown parts of the spectrum in addititon to better 
stability. 



It has been known since the early days of QCD that the mass splitting of hadrons 
in an isospin multiplet arises from two sources: virtual photon exchange, which 
can be related to virtual photon scattering by the Cottingham formula [1] and 
isospin symmetry breaking which arises from the difference between the masses 
of the light quarks m u and rrid as well as the difference between the light quark 
condensates (uu) and (dd). 

Let us start by a reminder of how electromagnetic and strong parts of the 
mass splittings arise [7] 

The quark mass term in the QCD Hamiltonian density reads 



The effect of turning on electromagnetism is described by the additional 
term 



2 The Calculation 



H{x) 



m u uu + rriddd + ... 



(1) 



e 2 f d 4 q f 




(2) 



which accounts for virtual photon exchange. 



is the electromagnetic current expressed in terms of quark fields. 
Following Cottingham[3] eq. (TTJ) can be transformed into 



(3) 
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H em (0) = ±^ J dq oy /Q*-qli J dtxe^Tj^x)^) 

= ^2 [ % j dq jQ*-<$.i J d'xe^Tj,(x)UO) (4) 
+ £2/°°^ f Q Q dq ^Q 2 -qli J d*xe^TjV(x)i p (0) 

where Q 2 — —q 2 and A is a divider between small and large values of Q 2 . 

Matrix elements of the low energy part of the low energy contribution to eq. 
([3]) can be obtained from experimental observation, the insertion of intermediate 
states, etc., this term we call H* oul (0), the "Coulomb" part. In the (divergent) 
high energy part we can use the operator product expansion (OPE) of the 
electromagnetic currents with the result 

e 2 A 2 - e 2 A 2 
flem(0) = H coul (0) + ^-5^—-mddd- — ^ In — .m u uu + ... (5) 

with A an ultraviolet cutoff. 

When this is added to expression (JTJ the divergent terms are absorbed in 
the renormalized masses and the Hamiltonian density takes the form 

H(Q) = H cou i(0) H .{uu + dd) H \uu - dd) (6) 

The first and last terms in the equation above break isospin symmetry and 
contribute to the neutron-proton mass difference SM np . 

ovn 

6M np = (p \H coul (0) I p) + — . ( P \u 3 1 p) = SM coul + 5M q (7) 

with 



2M P 



Srn = rrid — m u 



U3 = uu — dd (8) 

This separation depends of course on the value of the divider A but this 
dependence is extremely weak: the relative change in the value of SM q for two 
values of Ai and A2 is 

A 2 

- 1.5* KT 3 ln-f (9) 

5M cou i has been estimated a long time ago [5\. A recent evaluation [6] which 
is here adopted is 

8M coul = (-1.30 ± .50)AfeV r (10) 
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The measured value SM np — l.29MeV yields then 

{SM q ) cxp = (2.60 ± .50)MeV (11) 
The theoretical expression is given by eq. ([7]) 

SM q = 8mU (12) 

with 

U=(p\u 3 \p)/2M P (13) 

The aim of the present calculation is to evaluate U using the asymptotic 
forms of the various matrix elements given by QCD and the incomplete available 
knowledge about the nucleon continuum with an emphasis on the stability of 
the calculation. For this purpose I follow the external field approach of Ioffe 
and Smilga [5] in which the quarks are coupled to a weak external scalar field 
S(x) through an additional term to the QCD Lagrangian 

AL = —S(x).(uu - dd) (14) 

S(x) can be taken a constant. The correlation function of the nucleon cur- 
rents in the presence of S(x) is 

U(S,q) = t J d 4 xe lq - x (0\T Vp (x)r^(0)\Q) s (15) 
where rj p (x) is the proton interpolating field of Ioffe [2 J 

r] p (x) = eafc c [/J(x)C7 M l/fc(x)757 /J d c (a;) ( 16 ) 

where a, b, c stand for colour indices and C — — C T is the charge conjugation 
matrix. Lorentz covariance and parity allow the decomposition 

U(S,q) = U l (S,q) +7-«n«(S, q) (17) 
To first order in S the two invariant functions can be written as 

U i (S,q) = Ui(q 2 ) + SU\(q 2 ) (18) 

where IIq denote the invariant functions in the absence of the external field 
and U\ are the linear responses to the latter which, in QCD, can be expressed 
via the OPE to various vacuum condensates. 

The external field will contribute in two different ways: by directly cou- 
pling to the nucleon fields which enter in the nucleon current and by polarizing 
the QCD vacuum. This will introduce a susceptibility \ which describes the 
response of the quark condensates to S, 

(uu)s — (uu) — xS(uu) 

(dd} s = (dd) + x S(dd) (19) 
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where (qq) — (0 \qq\ 0) 
Using eq. (|13p one obtains 



X (uu) = -l dtx(0 \Tu 3 (x)u a (0)\ 0) (20) 
and, treating the term ^Smu^ in the QCD Lagrangian as a source term yields 

El 



X5m = -7 + 0{5m z )= - 7 (21) 

(uu) 

sates. 



with 7 = -^jj- — 1 a measure of isospin violation in the light quark conden- 



The QCD expression for Ilf (t) has been evaluated in [8] 

n?(*) = (coM-t) + ^ + ^ + ...) (22) 

with 
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Ci = -X« 2 )-^(W), (23) 
C2 = (x + Xm)mo(g§) 2 

X m denotes the response of the gluon condensate to the external field and, to 
a good approximation x-m = x|H]i 

m 2 = (9sq£^q)o ~ {M ± 2Q)Gey2 (24) 

At low energies II* (i) has double and single isolated poles at the nucleon 
mass squared as well as a cut on the positive real axis starting at t t h = (M p + 

-2XHl p .U A 

Ul{t) ~ (t-M*r + {t-M>) + ~ (25) 

where A is the coupling of the nucleon to its current 

(0|77 p |p) = XU P (26) 

In order to relate the residue U to the QCD condensates consider the integral 
of the product Ilf (t).F(i) over a contour C in the complex t plane. The contour 
C consists of two straight lines parallel to the cut, sandwiching it from above 
and below, and running from threshold to a value R and a circle of radius R 
taken large enough to allow the replacementof II' (t) by its asymptotic form eq. 
([H]) on it. 



F(t) is a so far arbitrary entire function. Using next Cauchy's theorem we 
obtain an expression relating the residue at the pole to an integral over the cut 
plus an integral over the circle where the QCD expression for the amplitude can 
be used. The integral over the cut cannot be calculated as the information on 
the integrand over the continuum is scarce. We shall use the arbitrariness in 
the choice of F(t) to minimize this contribution so it can be safely neglected. 

First take 

m = { l -il) m (27) 

in order to eliminate the contribution of the simple pole which represents 
nucleon to continuum transitions. Cauchy's theorem then yields 



2\ 2 M P U.P(M 2 ) = ~ j dtF(t)ImUl(t) + j) dtF{t){n\{t)) QCD (28) 

The second integral in the expression above runs over the circle of large 
radius R. 

The coupling A itselfcan be obtained from the nucleon mass sum rule PJ [3J [9J 
using P{t) as an integration kernel with the result 

- M P X 2 P(M 2 ) = - [ dtP{t)ImH 2 (t) + ^~S dtP(t)Il^ CD {t) (29) 

T Jth 2l7T J 

with 

Ur D (t)=B 3 tH-t) + ^ + ^- + ... (30) 



B 3 = i ^(^)(l+2-«) 
B 7 = -^(qq}(aGG) 

B 9 = Att 2 .lS6a(qqf (31) 

a = ag/n is the strong coupling constant. The ratio of eqs. (|28[) and (|29D 
yields the residue U of interest. 

Not enough experimental data is available to allow the evaluation of the first 
integrals on the r.h.s. of eqs. (|27p and (|28p. only the positions of the i + and \ 
resonances are known and the background is impossible to model realistically. 

The choice of the function P(t) aims at reducing this contribution as much 
as possible in order to allow its neglect. In the vast domain of QCD sum rules 
the usual choice would be P(t) = cxp(— t/M 2 ) where the magnitude of M 
(the Borel mass) determines the strength of the damping of the contribution of 



G 



the integrals over the continuum. If M is small the damping is good but the 
contribution of the unknown terms in the QCD asymptotic expansion of the 
amplitudes increases rapidly. If M increases the contribution of the unknown 
terms decreases but the damping worsens. An intermediate value of M has to 
be chosen from stability conditions which are not always met. 

I advocate here another choice for the integration kernel which stems from 
the observation of the spectrum of the positive and negative parity nucleon 
resonances N(1440), N(1535), N(1650) and N(1710). The position of these res- 
onances imply that the main contribution to the integrals over the nucleon 
continuum arises from the interval I : 2.0GeV 2 ~ < ~ S.OGeV 2 . The choice I 
adopt for the damping polynomial is 

P(t) = 1 - a[t - a' 2 t 2 (32) 

The coefficients a[ and a' 2 are chosen so as to practically annihilate the 
contribution of the continuum on the interval I, e.g. they can be chosen so as to 
minimize the integral L dtP(t) 2 . This choice reduces the relative contribution 
of the integrals over the continua to only a few percent of their initial values 
and allows their neglect. Numerically 



a[ = M7GeV- 2 , 

a' 2 = -.160GeV- 4 (33) 

We have thus 

F(t) = l-ait- a 2 t 2 - a 3 t 3 (34) 

with 



oi = 1.947GW" 2 , 
a 2 = -1.084GW" 4 , 

a 3 = .\84GeV- 6 (35) 

If the contribution of the nucleon continua are then neglected eqs. (|28p and 
(|2~g|) become 

2\ 2 ^-P(M 2 ) = C (l -°±R- |i? 2 - |i? 3 ) +C 1 - ai C 2 (36) 

X 2 M p P(M 2 ) = -B 3 R 2 (^ - yP- yi? 2 ) - B 7 + a[B 9 (37) 
the ratio of which yields U. 

The choice of R is dictated by stability considerations, it should be large 
enough to allow the use of the QCD expressions in the integrals over the circle 
of radius R but not too large to invalidate the neglect of the integral over 
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0.95 




R in GeV 



Figure 1: The right hand sides of (|36|) (A) and of (f37|) (B)as a function of R 
(not normalized) 



the continuum. The optimal value of R is chosen in the stability region of 
both equations (|3l)l) and (|37l) Indeed, as can be seen in Fig. 1 both terms are 
practically constant for 1.5GeV 2 ^ R ^ 3GeV 2 . 



3 Evaluation of the susceptibility x 

The next task is to calculate the susceptibility \ which enters in the expressions 
for the condensates. For this purpose I use a simple model independent method 

P33- 



Consider the correlators 

Mf(t = q 2 ) = i f d 4 xe^ x (0 TD l3 {x)D] t {Q) 0) 



(38) 



whereAj(a;) = (rrtj Trrii) : q 3 (x)i(^)<li(x) : 

is the divergence of the vector (axial- vector) current. At zero momentum a 



Ward identity relates \t2p(0) to the quark condensates [IT], viz: 



< 5) (0) = -(m 3 ± m^^f, T tf^n 



(39) 



we have therefore 
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-{m d - m u )(dd - uu) = * ud (Q) 

~(m d + m u )(dd + uu)o = y 5 ud (0) (40) 

Or x = -2(m d + m u )V(°)/* 5 (°) with #(i) = (m d - m u ) 2 ip(t). 

Consider first the scalar correlator which exhibits a cut on the positive t axis 
running from threshold to oo and a contourC similar to the one used previously 
for the nucleon correlator. The integral over C of the quantity i/>(t)/t gives i/j(0). 
As before the integration kernel is modified so as to minimize the contribution 
of the + continuum, i.e. consider the integral 



c 



dtA(t)ip(t), (41) 



A(t) = ~-b -b 1 .t (42) 

we have 

$(0) = - [ dtA(t)iJ>(t) + (t dtA(t)1>(t) (43) 

7T J cut 2«7T / 

In the first term on the r.h.s. of the equation above the integrand is impos- 
sible to model realistically and mostly unknown except for the existence of two 
isovector + resonances at m\ = .97GeV 2 and m\ — 2.lOGeV 2 . It is expected 
that the overwhelming contribution to this integral over the cut is provided by 
the interval .90Gev 2 ^ t £ 2.20GeV 2 

The choice of the constants bo and b\ such that A(mf) = A(m 2 ) — will 
practically annihilate the contribution of the integral over the continuum so that 
eq. (|42l) becomes 

^)-^-j> dt ^QCD{t)A{t) (44) 

As for the pseudoscalar correlator it can be treated in a similar fashion , the 
only difference being the presence of a pion pole which is by far dominant so 
that 

* 5 (0) - 2f 2 ml (45) 

Eqs. (gOl) to (g3D yield then 



W _ j = ( m i- m l ) r 

(dd) f 2 ml 



7=^-1= (46) 



with 

iP(q 2 ) = i f d 4 xe^ x (0 \T(q(x)q(x))(q(0)im\ 0) (47) 
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and give finally for the susceptibility 



K + m u ) 

The perturbative part of the correlator eq. (|4"7) is known to 5 loops in 
addition to 2 non-perturbative terms [T2"] 



il> QCD {t) = A t\n{-t) + ^- + ^ + ... (49) 



where 



A = -^(1 + ^-as + 14.1785<4 + 77.3535a 3 + 511.696a| + ...) 
A 1 = -l(a s GG)(l + ^. a (-lGeV 2 ) + ^o) 

^2 = ^n 2 a s (qq) 2 (50) 

So that eq. (gj) gives 

V(0) - A (i? - yi? 2 - -i? 3 ) - Ml - Ma (51) 

Once again the expression above is stable in the interval 1.5GeV 2 ^ R ^ 
2.5GeV 2 . The value of the susceptibility follows: 

X = 1.03GeF _1 (52) 

4 Results and Conclusions 

Before joining pieces together let us look at the values of the condensates which 
enter in the Ai, i^and C\ : 

(qq) is obtained from the Gell-Mann, Oakes, Renner relation[13j with the 
values of the quark masses m u = (2.9 ± .2)MeV,m d = (5.3 ± A)MeV [H] 

^ 2 {qq) = -(-79 ± .Ql)GeV z (53) 
The standard value is taken for the gluon condensate 

(aGG) = (.012 ± .006)GeV 4 (54) 

and 

((qq) 2 ) = n(qq) 2 (55) 

where k quantifies deviations from factorization. 
We have finally 
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A/f 1 ' ' 

U = _^(C i?(l-^i?-|i? 2 -|i? 3 )+C 1 -a 1 C 2 )/(-i?3i? 2 (i-|i?-^i? 2 )-i? 7 - a ' 1 i? 9 ) 

(56) 

Let us look at the sources of uncertainty in the expression above. 

The first stems from the incomplete knowledge of the OPE: with our choice 
of polynomials the numerator of eq. (|56[) should be augmented by two additional 
terms (I2C3 +036*4 and the denominator by a' 2 Bn. The higher order condensates 
C3, C4 and Snare of course unknown but it is possible to estimate them using 
the method of Pade' approximants. It turns out that the error introduced by 
the neglect of the higher order condensates amounts to no more than 2%. 

Another source of uncertainty stems from the choice of the coefficients a[ , a! 2 
(and consequently of the coefficients a\, ai-, 03 )., varying these within reasonable 
limit in order to deplete the contribution of the interval 2.0GeV 2 ~ t ~ S.OGeV 2 
introduces an uncertainty of ~ 6 — 7%. in the value of U. 

An additional uncertainty of course is the one arising from the values of the 
condensates themselves. 

All these uncertainties however are overwhelmed by the one coming from the 
undeterminacy in the value of k which measures the deviation from factorization 
in the value of the four quark condensate. There is no consensus on the value 
of this quantity in the litterature. Phenomenological estimates Q~5] place it in 
the range 1 ^ k ^ 4. 

With a central value [14] 6m = 2 AM eV, 5M q is plotted against k in Fig. 2 
where the experimental limits are also shown. 

Taking in consideration the uncertainty in 8m we see that the experimental 
value of SM q is reproduced for 

1.7<k<2.3 (57) 

It is instructive to compare the results obtained to the ones given by use of 
the exponential damping kernel P(t) — e ~aP". Expression ([56)) for U is then 
replaced by 

M 2 f M 2 M 2 rR/M 2 p 

U = -f{C M± J da;e -(x-^)-M p 2 C 1 + (l+^)C 2 )/(i?3M 4 y o dxxe-x+Br-jJL) 

(58) 

which reproduces qualitatively the same results as before but which fail to 
show good stability in variations of the Borel mass M 2 as can be seen from Fig. 
3 

To conclude, I have presented an evaluation of the neutron-proton mass 
difference SM np using polynomial kernels in dispersion integrals taylored to 
reduce the contributions of the unknown parts of the continua to an extent 
which allows their neglect and which guarantees the stability of the calculation. 

The numerical result depends sensitively on the value of Kwhich quantifies 
the deviation of the value of the four quark condensate from the one given by 
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Figure 2: The strong part of the neutron-proton mass difference 5M q as a 
function of k 
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the factorization assumption and reproduces the experimental value of 5M, 
forK ~ 2. 
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